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Abstract 

In the present article we liave found tlie complete energy spectrum and tlie 
corresponding eigenfunctions of tlie Dirac oscillator in two spatial dimensions. 
We show that the energy spectrum depends on the spin of the Dirac particle. 
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Recently, Moshinsky and Szczepaniak [|l| have proposed a new type of interaction in the 
Dirac equation which, besides the momentum, is also linear in the coordinates. They called 
the resulting the Dirac equation the Dirac oscillator because in the non-relativistic limit it 
becomes a harmonic oscillator with a very strong spin-orbit coupling term. Namely, the 
correction to the free Dirac equation 

^^ = (/97P + «^c (1) 

reads 

p ^ p — imuPr (2) 

after substituting (|]) into (|l|) we get an Hermitian operator linear in both p and r. Recently, 
the Dirac oscillator has been studied in spherical coordinates and its energy spectrum and the 
corresponding eigenfunctions have been obtained . A generalization of the one-dimensional 
version of the Dirac oscillator has been proposed by Dommguez-Adame. In this case, the 
modification of the free Dirac equation, written in cartesian coordinates, is made by means 
of the substitution m 'm — i'~f^'~f^V{xi). Obviously, for V{xi) = muJXi we have the standard 
Dirac oscillator. Here, as well as for the three dimensional Dirac oscillator 0, bound states 
are present. 

An interesting framework for discussing the Dirac oscillator is a 2+1 space-time. The 
absence of a third spatial coordinate permits a series of interesting physical and mathematical 
phenomena like fractional statistics and Chern-Simmons gauge fields among others. Since 
we are interested in studying the Dirac oscillator in a two-dimensional space, a suitable 
system of coordinates for writing the harmonic interaction are the polar p and coordinates. 
In this case the radial component of the modified linear momentum takes the form: pp — 
imuj(3p. It is the purpose of the present paper to analyze the solutions and the energy 
spectrum of the 2+1 Dirac oscillator expressed in polar coordinates. 

One begins by writing the Dirac equation (|1|) in a given representation of the gamma 
matrices. Since we are dealing with two component spinors it is convenient to introduce the 
following representation in terms of the Pauli matrices ||5| 
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= (Ji, /372 = 50-2, /3 = (T-i 



(3) 



where the parameter s takes the values ±1 (+1 for spin up and —1 for spin down). Then, 
the Dirac (0) equation written in polar coordinates reads 



with, 



a Op + a [ mLAjp) + ta m 



(4) 



where the spinor \E' is expressed in the (rotating) diagonal gauge, related to the cartesian 
(fixed) one by means of the transformation S{p, ??) [0] 



(5) 



where the matrix transformation S{p, 'd) can be written as follows 
noticing that S{p,^) satisfies the relation 



(6) 



S{p,^ + 2n) = -S{p,^) 



(7) 



we obtain, 



^(t9 + 27r) = -^(^9) 



(8) 



so we have = N + 1/2, where is an integer number 

Using the representation (^), the spinor equation (^ can be written as system of two 
first order coupled differential equations, 

d k^s 



i(E — m)'^i = (- — I imujp)'^2 

dp p 



(9) 



i{E + m)\E'2 



, d 
'dp 



Lis 



+ imujp)'^\ 



(10) 



where 

substituting (|^) into (||) and vice-versa we arrive at 
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T 1) 2 2 2 , 2 



(12) 



It is not difficult to see that the solution of of the second order equation (|T^) for \E'i can be 
expressed in terms of associated Laguerre polynomials L\.{x) [|^|[ as follows 

^1 = ci exp(-x/2)x5(i/2+^)L;^(x) (13) 

where we have made the change of variables 

X = mujp^ (14) 

where p satisfies the relation 

p = ±{hs - 1/2) (15) 
and the natural number n satisfies the relation 

_ ^2 

+ {lTl){2hs-l) = 4n (16) 



Since the function must be regular at the origin, we obtain that the sign of p in (|T3D is 
determined by the sign of s. In fact, for k^s > we have that \E'i reads 

^1 = ci exp(-x/2)x'=^^/2^^^.-i/2(^^) (^7) 
substituting (0) into (p!OD we arrive at 



v]>2 = 2^ci^^:^^^^ exp(-x/2)x('^^^+i)/2^^\+'/'(a;) (18) 
E + m 



where ci is an arbitrary constant. 

Analogously, we obtain that the regular solutions for k^s < are. 



^2 = C2 exp{-x/2)x-'''''/^L-''^'-^/\x) (19) 



E -\- m 



where ci is a normahzation constant The expression ([16|) can be rewritten as follows 

-m^ =4:[n- Q{~hs) {hs - 1/2)] mu (21) 

where B(x) is the Heaviside step function. Then from the relation (^) it is clear that the 
energy spectrum of the 2+1 Dirac oscillator depends on the value of s. Notice that for 
positive values of k^s there is not degeneration of the energy spectrum. For k^s < we 
observe that all the states with (n ± /, k^s — 1/2 ± /), where / is an integer number, have the 
same energy. In this direction there are some differences with the spherical Dirac oscillator 
0. Despite in both cases bound states are obtained, for the 2+1 Dirac oscillator the energy 
spectrum is degenerate only for negative values of k^s. In order to get a deeper understanding 
of the dependence of the energy spectrum on the spin we can take the nonrelativistic limit 
of the Dirac equation (^. In order to do that, it is advisable to work with Eq. (|T2|). The 
Galilean limit is obtained by setting E = m + e, and considering e << m. Taking into 
account that the first two terms in Eq.(|l^) are associated with the operator P^, we obtain 
in the nonrelativistic limit 

uj(k^s ± - + — = e 22) 

2m ^ 2^ 2 ^ ^ 

notice that Eq. corresponds to the Schrodinger Hamiltonian of a harmonic oscilla- 

tor with an additional spin dependent term given by —uj{k^s ± |). This contribution is 
proportional to the frequency of the oscillator. 

It would be interesting to analyze the Dirac oscillator in more complex configurations 
where electromagnetic and gravitational interactions are present, regretfully in this direc- 
tion the possibilities of finding exact solvable examples are limited to those were the Dirac 
equation with anomalous moment are soluble. A detailed discussion of this problem will 
be the objetive of a forthcoming publication 
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